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Abstract 

A remarkable extension of Rayleigh-Schrodinger perturbation method is found and 
described. Its (N + q) x (N + 1)—dimensional Hamiltonians are assumed emerging 
during quasi-exact constructions of bound states. At all q > 1, the role of the 
traditional single eigenvalue is taken over by an energy/coupling q— plet. In a way 
circumventing both the non-linearity and non-Hermiticity difficulties, the corrections 
are defined by compact, q —dimensional matrix-inversion formulae. 
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1 Introduction 


The existence and efficiency of several sufficiently simple perturbation techniques 
represents one of the not entirely negligible reasons why the abstract formalism of 
Quantum Theory finds so numerous and impressive practical applications. Usually, 
one is given a “realistic” (i.e., quite often, fairly complicated) Schrodinger equation 

(H-E) \iP) = Q 

and succeeds in its approximative replacement by a significantly simpler equation 

(i/ (0) - U (0) ) |^ (0) ) = 0 

where 

(1) the difference between the “correct” and “approximate” measurable quantities 
(here, between energies E n and E at a few level-indices n) proves small, 

(2) the latter, “zero-order” problem is comparatively easy to solve, and 

(3) a suitable (e.g., Rayleigh-Schrodinger [1]) version of the so called perturbation 
theory is found as a source of a recipe for a systematic improvement of the 
approximation. 

In a broader perspective, one may keep the energy E fixed and switch to a modified 
Schrodinger operator, H — E = H\E ) — gV = H" — S(E, g) = . ... An appropriately 
modified perturbation theory might then serve as a tool for the determination of the 
“fine-tuned” couplings g [2] or, whenever necessary, of the whole multiplets of the 
couplings <?i, g 2 , ■ ■., g q -1 which may even re-incorporate the energy itself by putting, 
formally, E = — g 0 . In the latter setting the “original” or “realistic” problem quite 
often acquires the “q— plet-eigenvalue” generalized form 

[H - S(g 0 ,g 1 ,...,g q _ 1 )} \g 0 , g u ..., g q -i) =0 (1) 

which will be considered and studied in what follows. 


2 Polynomial-oscillator sample of equation (1) 

Differential Schrodinger equation 

+ V(x) tp(x) = Eip(x) (2) 
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for a particle moving in a one-dimensional polynomial well 


V(x) = giX 2 + g 2 x 4 + ... + g 2q+ ix 4q+2 , g 2q+ i > 0 (3) 


may be used as an illustrative example. Its solutions prove often tractable via the 
infinite Taylor series 

OO 

if(x) = exp[— P(x)] x y h n x 2n+p (4) 

»i=0 


where the integer p = 0 or 1 characterizes the parity. It may be re-interpreted as an 
angular momentum t = p — 1 = 0,1,... in three dimensions, etc. 

In place of the general exponent P(x) we shall only use here a polynomial 

2fc+2 


P (x) = E 


X 


k =0 


2k+ 2 


fk ■ 


( 5 ) 


The asymptotically optimal choice of its coefficients fk is uniquely determined by 
the WKB condition 

V(x) = [P'{x)} 2 + O (x 2q ) (6) 

or, if you wish, 


<? 

92q+l = ffi 92q = ^fqfq—lj ■ ■ ■ i 9q+ 1 = fq~jfj (^) 

3=0 


which hxes all the q +1 coefficients fk as functions of the first g+1 dominant couplings 
92 q+i, 92 qi ■ ■ ■ , 9q +1 in potential (3). Thus, we have /o = ^[ffi for harmonic oscillator 
(i q = 0) or /i = and /o = g 2 /(2fi) for the sextic anharmonic oscillator (q = 1) 
etc. 

The insertion of (4) transforms our differential Schrodinger eq. (2) + (3) into the 
infinite-dimensional matrix problem (1) where the non-square matrix of the system 
possesses the (q + 2)—diagonal banded form, 


H ~ S(g 0 , <7i,... ,g q -i) 


Bo 

Co 



4 1 ’ 

B x 

c, 


•• 


olW B, 

Cq 



A { °1 ■■■ 

B m C m 


( 8 ) 
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The g/j—dependence of the matrix elements is transparent, 


C n — (2 n + 2) (2 n + 2p + 1), B n — —go — fo (4n + 2p + 1) = A®, 

A^ = -/i (4™ + 2p - 1) + / 0 2 - ^1, ^i 2) = -f 2 (4 n + 2p-3) + 2/ 0 /i - # 2 , 

(9) 

Ai^ = — fq (4w + 2p + 1 — 2g) + (fofq-l + flfq-2 + • • • + fq-lfo ) — fl'q, 

n = 0 , 1, ... 


(remember that — g 0 = E is the energy). We have 

q 

S(go,gi,...,g q ~i) = 9z-^i ( 10 ) 

£=1 

where the generalized vectorial eigenvalue g = {go, gi, ■ ■ ■, g q -i} is attached to a 
q —plot of auxiliary one-diagonal unit-like infinite-dimensional matrices 


2.1 Illustration: q = 0 (harmonic oscillator) 

At q = 0 we may re-scale the coordinate x and set V (. x) = x 2 and g± = fo = 1 giving 
ip{x) = exp(— x 2 /2) x ... etc. From the formal secular equation 


det [H - S(g 0 ) 


det 


( Bq Co 
0 B\ Ci 




= 0, q = 0 


J 


we deduce that the harmonic-oscillator spectrum is reproduced once we let the secular 
determinant vanish by setting one of its diagonal elements B m equal to zero. Thus, 
the q — 0 example may be skipped as trivial in what follows. 


3 Magyari’s finite-dimensional equations (1) 

In a pioneering letter [3], E. Magyari discovered the existence of the so called quasi- 
exact [4] bound states in some special cases of the general polynomial potential 
wells (3). In our present language, his claim is that the infinite-dimensional matrix- 
equation representation (1) + (8) of the polynomial oscillators admits a rigorous 
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and exact finite-dimensional truncation. We only have to select the very specific 
truncation-compatible “intermediate” coupling 


9q ~ ~fq (4-/V + 2q + 2p + 1) + (fofq-1 + flfq-2 + • • • + fq-ifo) (11) 


in order to arrive at the finite-dimensional matrices in eq. (1), 


H — S — 


( B 0 C 0 
4 1} B 1 


A {q) 

q 


Cl 


A? B, C, 



aO /td) p 

/i N-2 ‘ ‘ ‘ /1 V-2 n N- 2 

1 

to 



A (q) 

Wv-i 

-Bat-1 

Cn-i 




Bn 

V 

A fa) 

-rt-N+q-2 

-rt-N+q-2 

4(9) 

^N+q-l 

^N+q-l ) 


( 12 ) 


At every choice of the integer N = 0,1,2,... it makes sense to contemplate the 
next, O (x 2q ) order of precision in our WKB formula (6). This leads to eq. (11) 
and connects the value of the subdominant coupling g q with the x 2N+p —dominated 
truncated wave functions (4). 

We get A$ = 4 f q (N + q — n) and reduce the original infinite-dimensional linear 
algebraic system (1) to its finite-dimensional subset which consists of the mere N + q 
rows. As long as we assumed that hv+i = hjq +2 = ... = 0, only the IV + 1 Taylor 
coefficients h 0l hi ,..., hjy remain unknown and must be determined by eq. (1). In 
this sense the latter equation is overcomplete and, in effect, non-linear. The operator 
H — S becomes represented by a non-square, (N + q) x (N + 1)—dimensional matrix. 


3.1 Sextic quasi-exact oscillators (q = 1 and all N = 0,1,...) 

At the first nontrivial q = 1 we may re-scale the coordinate x and pick up, say, the 
even parity p = 0. This gives V ( x ) = x 6 + 2f 0 x 4 + (/q — 3 — 4N)x 2 and the bound 
states ^(x) = exp(— x A /4 — fox 2 / 2) x .... One arrives at the square-matrix secular 
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equation of the N + 1—dimensional form 


l B 0 Cq 

4 1} B 1 C ! 


det [H — *S'(^o)] — det 


B N -i Cn-i 


A 


(i) 

N 


= 0 , 


9 = 1 - 


B n ) 


The standard square-matrix techniques may be employed so that also all the quasi- 
exact constructions of the sextic oscillators with q — 1 may be skipped here as a 
mere routine exercise. 


3.2 Explicit solution at N = 0 and all q 

Once we contemplate the higher degrees q — 2, 3,... in our quasi-exact polynomial- 
oscillator Schrodinger-Magyari non-square-matrix bound-state problem (1) + (12), 
we must admit that their solution represents a challenge. The only manifestly non- 
numerical example occurs at N = 0 where the solution remains elementary, 


E = -go = /o(2p+ 1), 0i =/o-/i(2p + 3), 

fib-1 = (fofq-2 + flfq—3 + • • • + /q- 2 / 0 ) — fq- 1 (2 p +2 q — 3), N = 0. 

The asymptotically dominant couplings s^+i, g 2 q , • • •, g q +1 in the potential (3) remain 
freely variable while the q + 1 elements of the remaining set g 0 , g±,..., g q become 
fixed as their explicit quadratic functions. One may preserve such a pattern in all 
the N > 1 cases. 

3.3 A toy decadic example: q = 2 and N = 2. 

Before we get to the core of our present message and formulate the general pertur¬ 
bation method of solving eq. (1) at the arbitrarily large integers q > 2 and iV > 1, 
let us emphasize that the idea of the method may be most easily illustrated by its 
application to the elementary decadic models with q — 2 and with the re-scaled x 
such that /2 = yjgi = 1. The two free parameters f\ (= fiq/2) and /o [= (r/3 — ff)/ 2] 
define g 2 = 2/ 0 /i — 2p — 13 [cf. eq. (11)] and exponent P(x) [cf. eq. (5)]. They also 
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occur in the matrix elements of eq. (1) which reads 


1 -go - fo - 2p/o 2 + 4 p 0 N 

fo ~ .9i - 3/i - 2pfi -g 0 - 5/ 0 - 2 pf 0 12 + 8p 

8 fo ~ .9i - 7/, - 2p/, —fi'o — 9/ 0 — 2p/ 0 

V 0 4 fo - 91 - n/i - 2p/i / 


(13) 

and which should determine the two unknown eigenvalues go (representing the nega¬ 
tive energy) and (+ (= the non-variable coupling at a; 2 ) as well as the two coefficients 
ho and hi in the wave function (note that h? = 1 is a mere fixed normalization 
constant). 


4 Perturbation forms of Magyari’s equations (1) 

Studies of analytic potentials found a mathematical encouragement in the Regge’s 
complex angular momentum method [5]. In this context, our parity parameter could 
be prolonged to the complex angular momentum £ = p — 1 G C. Rather unexpectedly, 
the resulting extension of the range of p will enable us to simplify the bound-state 
problem (1) + (12) in the limit of the large \p\ S> 1 [6]. Let us now outline the two 
alternative possibilities of such a simplification in some detail. 


4.1 An elementary large —£ model 


In the toy example of paragraph 3.3 let us introduce the quantity A = l/(2p) and 
assume that the numerical value of this A may be chosen arbitrarily small. This 
implies that the identical decomposition 


H 


2 p 


H {0) + A LT (1) ' 


q = N = 2 


(14) 


of our Hamiltonian in eq. (13) with 


H W 


f ~fo 2 0 \ 

-/i ~fo 4 
0 -A -fo 

\o 0 -/, 


and H (l> 


f ~fo 2 0 \ 

/ 2 - 3/, —5/ 0 12 

8 / 0 2 - 7/, —9/ 0 

V 0 4 / 0 2 -ll/J 


acquires the manifestly perturbative character. The dominant H ^ is a very simple 
tridiagonal matrix while the more complicated and four-diagonal matrix remains 
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tractable as a mere perturbation. Thus, the current conditions of the applicability 
of perturbation algorithms are satisfied. We may expect that at the larger integers q 
and/or N, the same pattern will enable us to postulate the formal split of any given 
non-square Hamiltonian H in some formal Taylor series 

H = H (0) + A H (1) + A 2 H (2) + ... (15) 

where the choice of the “sufficiently small” parameter A is arbitrary. Even in the 
present letter this quantity may be re-defined as A = 1/(2 p + const), in the spirit of 
the popular “shifted large —£” expansions [7]. 

4.2 A more advanced large —£ arrangement of the same model 


T 


5 


/ h 0 \ 
S" 1 h X 

\ h '2 J 


= 0 


Let us abbreviate —go — 2pf 0 = 2pE, fo — gi — 2pfi = 2pF and re-scale eq.(13), 

( 2pE-f 0 2 + 4 p 0 \ 

2 pF - 3/i 2 pE - 5/ 0 12 + 8 p 

8 2pF — 7 fi 2pE — 9/ 0 
V 0 4 2pF — Ilfi ) 

with the “small” a = p~ l F in 

( ( 4 P) ^ 

u- 2 /(4p) 

a- 3 /(4p) 

u- 4 /(4p) ) 

This leads to another re-formulation of our problem (16), 

/ ho \ 

[H — S(s, t)} \s, t) = 0, j s,t) = S^ 1 hi 

\ h'2 ) 


(16) 


T = 


V 


and S = 


( 1 


V 


\ 


a 


) 


E F 

S —Ta’ *=~^ < 17 > 


and to another perturbation arrangement (15) of the Hamiltonian where the new 
unpert urb ed Hamilt onian 


H^ = 


( 0 

1 

0 \ 

0 

0 

2 

2 

0 

0 


1 

o ) 


(18) 
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is accompanied by the three non-vanishing one-diagonal matrix correction coefficients 
in eq. (15) with A = a — p -1 / 4 -C 1, 


/jl 1 ) — 

/ 0 0 0 \ 

3 0 0 

#(2) = _h 

/ 1 0 0 \ 

0 5 0 

, H i3) = - 

/ 0 i 0 \ 

0 0 6 

4 

0 7 0 

4 

0 0 9 

’ 2 

0 0 0 


V 0 0 11 ) 


V 0 0 0 ) 


V 0 0 0 ) 


The eigenvector | s,t) as well as the eigenvalue matrix S in eq. (17) should be also 
represented by the infinite Taylor-series ansatzs 


S = S^ 0) + A S (1) + A 2 S {2) + ... (19) 

and 

| s,t) = |s (0) ,f (0) > + A|s (1) ,f (1) ) + ... . (20) 

The key specific merit of such a sophisticated innovation lies in an unexpected and 
very significant simplification of the zero-order (i.e., \p\ —> oo limit of) Schrodinger- 
Magyari equations. For more details we may either refer to our older paper [8] or to 
the Appendix A where a few more technical details are summarized. These details 
are instructive since our particular q = N = 2 zero-order concrete model (18) admits 
a straightforward generalization to all N, at q < 5 at least [6]. 

5 The evaluation of perturbation corrections 

The generic infinite-dimensional matrix equation (1) may be solved, say, by the tech¬ 
nique of the so called vectorial continued fractions [9]. Perturbation re-arrangements 
of such a treatment of our problem exist but still carry a more or less purely numerical 
character [10]. Here, an entirely different approach is being employed. 

For the sake of definiteness, we shall only pay attention to the Magyari’s finite¬ 
dimensional constructions at arbitrary truncation integers N < oo. The scope of 
the method might be broader in principle but we decided to consider the Magyari’s 
problem only. Indeed, due to the practical difficulties with its solution it hardly 
possesses any applications beyond the first few qs or Ns at present. 

As we already indicated, our present approach to the Magyari’s problem (1) (we 
may call it “perturbed” problem) will be perturbative. We shall postulate 

\g 0 , gi,..., g q -i) = |4°\ 4°\ • • •, g q l i) + A | go\gi\ ■■■, 4- 1 ) + • • • ( 21 ) 

8 



and employ the parallel expansion (20) of the q —plet of eigenvalues which contains 
S = S{go, gi ,..., g q -i) and S {fe) = S(g 0 (k) , gi (k) ,gfl i) with k = 0,1,.... 

We may re-write our “perturbed” Schrodinger equation (1) order-by-order in A. 
One expects that for the study of similar equations, all the solutions are at our 
disposal not only for the overcomplete zero-order equation 


h { 0) \g ( o\gi\■ ■ -,A) = S(go (0) ,gi (0) ,■ • ■, A)l£o 0) ,3i 0) ,• • ■ ,s£-i) 


but also for its under-complete “left-action” partner 


) _(°) I rr(0) _ / / (0) (0) (0) i q/ (0) (0) (0) \ 

) • • ■ > 9q— 11 \£\9 0 : g 1 1 • • • i 9q— 11 h \9 0 i 9 1 i • • • i 9q— 1 / • 


On the level of the first-order corrections we complement the zero-order homogeneous 
equation (22) by the non-homogeneous linear algebraic problem using an abbreviation 
I 0) = |So° , .Si 0, .---.4-i>. 

'H m - S (0) ] IjS 11 , j! 11 , .... 4-i> = [s (1) - ff a) ] I 0). (24) 

Similarly we add its second-order descendant 

tf(°) - S (0) ] | g ( o\gf\ ..., 4- 1 ) = [S {1) ~ H (1) ] | 1) + [^ (2) - H^] \ 0) (25) 

(with abbreviated | 1) = > flg-i)) an( l so on - All the equations in this 

hierarchy have the same recurrent structure, 

'H (0) - S (0) ] | k) = | known (fc - 1} ) + S {k) \ 0) (26) 

with “input” |known (0) ) = -H {1) \ 0) at k = 1, etc. 

This completes our description of the overall strategy and algorithm. In depen¬ 
dence on the particular features of and of its eigenstates, a few more detailed 
technical aspects of our innovated recipe may be added in a more explicit form which 
is dependent of course. 

For our concrete construction of section 4.2 (which admits various generalizations 
[6]), the technical details are particularly rich and relevant and they contribute to a 
perceivable enhancement of the efficiency of the practical calculations. All of them 
may be found summarized in Apenclix B which offers a slightly more formal reading. 
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6 Discussion and summary 

6.1 The next-to-elementary quasi-exact bound states 

In various applications of Quantum Mechanics people sometimes forget that the 
occurrence of the compact and elementary “user friendly” wave functions is not 
restricted to the mere linear harmonic oscillator V{x) = x 2 . Due to Magyari [3] we 
know that the similar exceptional bound states, say, of the form 

^{qe){x) = x const exp(—polynomial (x)) x polynomial (x) (27) 

are much more generic and may be generated by virtually any polynomial potential 
at certain particular subsets of couplings and energies. 

Unfortunately, all his formally flawless generalizations of the q = 0 harmonic 
oscillator to q — 1,2,. .. found applications solely at q — 1 (cf., e.g., [11]). Certain 
q = 2 models attracted attention just at the very small N [12], In all the other cases, 
the unfortunate coincidence of the nonlinearity and non-Hermiticity of H proved 
deterring. In this sense, our present new approach might re-vitalize interest in the 
undeniable phenomenological merits of the quasi-exact models with the larger qs and 
Ns. 

In this direction we already re-analyzed a few particular non-square Magyari’s 
equations (1) at the freely variable dimensions N and discovered, purely empirically, 
that some of these equations exhibit certain unexpected “user-friendly” features up to 
q = 5 at least [6]. Moreover, it has been established by several groups of authors that 
in a way paralleling the exceptional q = 1 oscillators, some of their “first nontrivial” 
q = 2 descendants also admit the existence of quasi-exact multiplets which may 
be arbitrarily large [13]. Efficiently, their non-Hermiticity may be tamed by an 
appropriate modification of the “physical” inner product in Hilbert space [14], 

All these observations might revitalize interest of people in the Magyari’s non- 
square matrix problem (1). In our present paper we demonstrated that in spite of its 
non-linearity, rather surprisingly, this anomalous form of the Schrodinger equation 
still admits a systematic perturbative solution. We believe that the resulting facil¬ 
itated perturbative tractability might make Magyari’s oscillators significantly more 
attractive in applications. 
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6.2 Innovations in perturbation theory 

Perturbation analysis of nonlinear systems usually requires a narrow specification 
of their nonlinearity (cf., e.g., [15] for illustration). Hence, it was a nice feeling to 
reveal that the nonlinearity which characterizes the Magyari’s equations (1) still does 
not prevent us from employing the ideas of the Rayleigh-Schrodinger perturbation 
construction of their solutions. 

We should remind the readers that the use of the perturbation series is one of the 
most natural strategies, able to provide closed formulae as well as the upper and/or 
lower bounds of the measurable quantities [16]. Equally well it seems suitable for the 
study of the real and complex eigenvalues [17]. This puts our present perturbative 
approach to Magyari’s problem into a broader context and comparison. 

First of all, the specific merit and recommendation of our present approach for 
future calculations may be seen in the rapid growth of the complexity of the various 
non-perturbative alternative methods with the growth of N. Indeed, just the first 
few smallest integers N seem to have been considered in the “minimally nontrivial” 
q — 2 context up to now [18]. 

Secondly, our present method looks comparatively friendly also with respect to 
the growth of the integer q which measures the flexibility of the shape of the potential 
(3). In a way complementing and extending our thorough perturbation study at q — 1 
[19], we showed here that an apparent incompatibility of the perturbation-expansion 
strategy with the Magyari’s equations at the larger q proves false. We imagined that 
even some quite strongly non-square matrix forms of the Magyari’s Hamiltonian 
matrices with q > 2 remain tractable when both their left and right action are taken 
into account simultaneously. 

In this sense, a certain methodical gap has been filled by our present paper. 
Still, several open questions survive. For example, the finite-dimensional character 
of our present Hamiltonian matrices H might help to suppress the weight of some 
problems with convergence which often mar the use of the perturbation techniques for 
many infinite-dimensional H s [20]. In particular, one might recollect the encouraging 
rigorous “acceleration of convergence” results as obtained at the “trivial” q = 1 in 
ref. [19]. 
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Appendix A: Zero-order example (18) 


For the time being let us ignore the notation subtleties and changes of notation of 
section 4.2 and assume simply that a given toy unperturbed non-square Hamiltonian 
(18) may be assigned the generalized eigenvalues 


S < 0 ) (fff,g lb 


and eigenvectors 


( 9o <0 > 

9. (0) Jo (0) 


J. (0) Jo (0) 

9i (0) / 


= go 




I = 


( hf 
hf 

\ hf ) 


The exhaustive analysis of this problem would require a lot of space. Fortunately, it 
has already been performed elsewhere [8], with the most relevant result being that 
one just gets not more than two real solutions 


|M> = 

{ 1 1 
l 

1 

to 

1 

( 1 \ 
-2 


l i ) 

! 1 

l 1 / 


(28) 


They may be interpreted, as we already mentioned, as a doublet of the decadic 
quasi-exact large —i bound states. In some applications, one also needs to find all 
the complex eigenvalues [21] . Although they would still be tractable by our present 
perturbation recipe, we shall ignore these solutions here as unphysical. 


A.l. Left zero-order eigenvectors 

The left eigenvectors of our perturbed as well as unperturbed Hamiltonians H may 
be found of interest as solutions compatible with the conjugate equation 

9i, ■ ■ ■, 9q- 1 | H = (^(go, <7i, • ■ ■, g q - 1 | 9i, ■ ■ ■, 9 q - i) • (29) 

In a double-bra Dirac-like notation, these row vectors {^(go, 9i, ■ ■ ■ , 9 q -i\ are num¬ 
bered by an additional subscript £ = 1,2,... ,£, max . Such a convention emphasizes 
the ambiguity of the solutions. 
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Our special illustration hfl°) of eq. (18) offers the respective symmetric and 
antisymmetric left eigenvector 

(r(l, 1| = (1,1,1,1), <2<1,1| = (3, -1,1, -3) (30) 

at g 0 = 9i = 1, or 

(i(—2, —2| = (2, -1, -1, 2), ( 2 (-2, —2| = (0,1, -1, 0) (31) 

at r/ 0 — f)\ — —2. Arbitrary superpositions of these doublets may be employed in¬ 
stead. This reflects the incomplete (or rather “under-complete”) character of eq. (29). 


A.2. Reduced zero-order eigenvectors 

At any q > 1, one might prefer the choice of the row eigenvectors with as many zeros 
as possible. It is obvious that the linearly independent q —plets of the vectors of this 
type may be written as products 


{d9o ] , gi°\ - ■ ■, 9q~i\ = {j{e(g< 


( 0 ) 


,( 0 ) 


n (0) 
9 i ) 


n (0) ) 
• 9q-l) 


x% j — 1,2,... ,q (32) 


with (. N + 1) x (. N + q )—dimensional “reduction” matrices II j (= transpositions of 
some of the “unit-like” matrices with ^ j ^ in general). 

In our example at q — N — 2, the sample of vectors in eq. (30) may be replaced by 
the mutually conjugate “reducible” eigenvectors ( 3 ( 1 , 1 ! = (3,1,2,0) and (4(1, 1 1 = 
(0, 2,1, 3). Their explicit reductions 



( 1 

0 

0 

o N 


( 0 

1 

0 

0 

3 (1,1| = (3,1,2) x 

0 

1 

0 

0 

, ( 4 (1,1| = (2,1,3) x 

0 

0 

1 

0 


lo 

0 

1 

0 ) 


lo 

0 

0 

1 J 


may be abbreviated as (3(1, 11 = (i(p(l,l)|IIi and (4(1,1! = (2(^(1, 1 )| n 2 , respec¬ 
tively, with II j = Jj . 


A.3. Linearly independent sets 

The reduced, (N + 1)—dimensional “p—vectors” (j(p| may be interpreted as the left 
eigenvectors of the (IV + 1) x (N + 1)—dimensional square matrices [Ilj ( H 0 — S')]. 
These reduced eigenvectors must exist and be nontrivial since all their “secular” 
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determinants vanish. In such a context, one may sometimes need to guarantee the 
linear independence of the row vectors at the first q subscripts j = £ = 1,2 ,,q 
(say, by an appropriate re-numbering of their “overcomplete” available family). 

Sometimes, we may have to make the choice of IL,- = Jj, with j ^ j' for this pur¬ 
pose. For example, the mutually conjugate superpositions ( 3 (—2, —2| = ( 1 , 0 , — 1 , 1 ) 
and ( 4 (—2, —2| = ( 1 , — 1 , 0 , 1 ) of the states ( 31 ) lead to the respective factorizations 



( 1 

0 

0 

o N 


/ 1 

0 

0 

0 ^ 

I = (1, —1,1) • 

0 

0 

1 

0 

, ( 4 (— 2 , —2 = (1, —1,1) • 

0 

1 

0 

0 


lo 

0 

0 

1 ) 


\0 

0 

0 

1 J 


with the same form of the reduced array (i(q(— 2, —2)|. Such a degeneracy reflects 
merely the presence of the random zeros in the antisymmetric vector in (31). It 
may be avoided easily since the relation ( 2 (—2, —2| = (0,1,—1)11! = (1, —1,0) |n 2 
gives the other two alternative options for the second and safely linearly independent 
reduced vector ( 2 (^(—2, —2)|. 


Appendix B. A few technical aspects of our inno¬ 
vated perturbation recipe 

B.l. The explicit evaluation of the corrections to the per¬ 
turbed energies and couplings 

As long as our unperturbed Hamiltonian occurs also in eq. (23) which may be 
re-written in the form 

Gte(sS 0) , a[°\ • • •, gi°2 1 )I n, (i/ (0) - s (0) ) = o, j = 1 , 2 ,..., q , ( 33 ) 

we may multiply eq. (26) by the set of matrices n ,- and by the independent vectors 
(j{q (go°\ gi°\ ■ ■ ■, <4-i) I = (j(&l from the left. This leads to the system of relations 

0(^1 n i S(k) I °) = “0(^1 n i I known (fc_1) ) , j = 1, 2,..., q , (34) 

i.e., to an elementary q— dimensional matrix-inversion definition of the k— th-order 
energies/couplings (/^-i , 

E d,?"”’ 4-. = 4 1 ”' > i= 1,2..... 9- (35) 

f=i 
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Here, we merely introduced abbreviations <d kno " n) — — | known^ -1 ^) and 

^-.jknown) _ j-j-jjo) for scalar products of two known vectors in TV + 1 di¬ 
mensions. Their calculation may be facilitated by the evaluation of all the elements 

N+q 

[Gjtlmn = E [HiLk Wkn m , n = 0, 1, . . . , TV 

k =1 

in a preparatory step. One should remember that all the elements of both the non¬ 
square and sparse matrix factors are here equal to 0 or 1. The matrix elements of 
the product will often vanish, therefore. In the generic case with n., = Jj we 

even get 

(\j-Z\>N + l) =► jfJi = 0. (36) 

This means that the sparse matrix T of the coupling-defining system (35) has a 
band-matrix structure and its inversion is easier. 

B.2. Perturbed eigenvectors 

Whenever the solution of the zero-order eq. (22) is unique, the ambiguity of the 
vector | k) specified by the k— th eq. (26) lies in an uncontrolled admixture of the 
(TV + 1)—dimensional zero-order vector | 0). Once we construct any (normalized) 
basis in the corresponding (TV + 1)—dimensional space, 

(I /uiiL, 

and assume that the above vector coincides, say, with its zeroth element, | 0) = | /? 0 ), 
we may simply follow the textbooks and define the “right” projector 

N 

Qr — E l ft™) l 

m= 1 

which mediates the standard Rayleigh-Schrodinger normalization of | k) — Qr\ k). 

In parallel, we have to recollect that we already employed q independent rows of 
the k —th eq. (26) for the specification of via eq. (35). More specifically, we 
multiplied eq. (26) by q vectors (j(^| n^- from the left. As long as these vectors were 
chosen as linearly independent, they span a q —dimensional subspace in the “larger”, 
(TV + q )—dimensional vector space. This allows us to assume that the orthogonal 
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complement of this subspace is spanned by an N— plet of some (N + q)— dimensional 
vectors | a n ). These vectors define our second, “left” projector 

N 

Ql ^ ) I ®n) (p^n | • 

n= 1 

We are now ready to replace our non-homogeneous algebraic eq. (26) by its subsys¬ 
tem 

Ql [. H (0) - S (0) ] Q r | k) = Q l \ known^ -1 '*) + Q L S « | 0) (37) 

where the second term on the right-hand side has been made “known” in the previous 
subsection. By construction, the left-hand-side N by N matrix is invertible and we 
have 

I k) = Qr q l [ff(o) ]_ g(o)] q r Ql (| known (fc - 1} ) + S {k) \ 0)) (38) 

which is our final and compact “generalized Rayleigh-Schrodinger” explicit formula 
for the k— th correction to the wave function. 

B.3. Left perturbed eigenvectors 

We have seen that in the Magyari-equation context, the left vectors only played a 
role in their unperturbed , zero-order form so that the description of their pertur¬ 
bation construction might simply be skipped as redundant. Nevertheless, for the 
sake of completeness, let us add a few remarks also on the q —plct of the row (left) 
eigenvectors of the non-square Hamiltonians and on all the space spanned by their 
superpositions 

((^1 = X) Q • (z(9o,gi, ■ ■ ■,g q -i\■ (39) 

e=i 

By assumption, they all satisfy the perturbed conjugate eq. (29) and we may expand 
both their separate components and the coefficients in the power series in A. The 
combination of these series 

Q(9o, 9i, ■ ■ ■, 9g-i\ = Q(9o ] , 9i ] , ■ ■ ■, 9q-i\ + ^ Q(9o\ 9i ] , ■ ■ ■, 9q-i\ + ••• (40) 

and = J2k A k C'j: k) gives 

m = E A fc ((^% (W {k) \ = EE ct n) • I ( 41 ) 

k g=l n=0 
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where, in a naturally simplified notation, (n | = 

Let us now enter the most characteristic postulate of all the Rayleigh-Schrodinger- 
type perturbation recipes which removes the ambiguity of the eigenvectors by the 
consequent requirement of a complete absence of any zero-order component (^(0 | in 
any correction ((' ip( k l\ with k > 1. 

In our present case this option leads to the two independent consequences. Firstly, 
in the light of eq. (41) we have to put = C^ = . . . = 0 [i.e., n — k in eq. (41)] 
and get just C% = Cj?\ Secondly, all the vectors (^(n | must be constructed as 
perpendicular to all (j( 0 | at any n > 1, i.e., in our notation of section 6.2, 

(t( n \ = (t( n \QL, n — 1,2,, £ = 1,2 ,...,q. (42) 


With this normalization we may now insert eq. (41) in eq. (23) and arrive at an 
analogue of eq. (26). Of course, we may skip the repetition of the reconstruction of 
S ( fc ) and jump immediately to the following k > 1 analogue of eq. (37), 

(t(h | Ql [H W - S^} Q r = (z{ known^- 1 ) \Q R , £ = 1, 2,..., q . (43) 


With 

(*< known (0) | = ( ? (0 | [S {1) - H (1) ] (44) 

at k — 1 etc, the trivial N by N inversion gives the final formula 

1 

{n( k I = (O known | • Qr ^ ^ (0) _ g(0) j ^ Ql , £ = l,2,...,g. (45) 

The detailed account of our Magyari-inspired recurrent recipe is completed. 
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